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ABSTRACT 

We show how accretion rate governs the physical properties of a sample of unobscured broad-line, 
narrow- line, and lineless active galactic nuclei (AGNs). We avoid the systematic errors plaguing 
previous studies of AGN accretion rate by using accurate accretion luminosities {Lint) from well- 
sampled multiwavelength SEDs from the Cosmic Evolution Survey (COSMOS), and accurate black 
hole masses derived from virial scaling relations (for broad-line ACNs) or host-AGN relations (for 
narrow- line and lineless AGNs). In general, broad emission lines are present only at the highest accre- 
tion rates {Lint/ LEdd > 10"^), and these rapidly accreting AGNs are observed as broad- line AGNs 
or possibly as obscured narrow-line AGNs. Narrow-line and lineless AGNs at lower specific accretion 
rates {Lint/ LEdd < 10~^) are unobscured and yet lack a broad line region. The disappearance of the 
broad emission lines is caused by an expanding radiatively inefficient accretion flow (RIAF) at the 
inner radius of the accretion disk. The presence of the RIAF also drives Lint/ LEdd < lO^" narrow-line 
and lineless AGNs to 10 times higher ratios of radio to optical/UV emission than Lint/ LEdd > 10~^ 
broad-line AGNs, since the unbound nature of the RIAF means it is easier to form a radio outflow. 
The IR torus signature also tends to become weaker or disappear from Lint/ LEdd < 10"'^ AGNs, 
although there may be additional mid-IR synchrotron emission associated with the RIAF. Together 
these results suggest that specific accretion rate is an important physical "axis" of AGN unification, 
described by a simple model. 

Subject headings: galaxies: active — galaxies: nuclei — quasars: emission lines — accretion, accretion 
disks 
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Supermassive black holes (SMBHs) are now known 
to be ubiquitous in th e centers of all massive galaxies 
(jMagorrian et all 119981 ). SMBHs grow in an "active" 
phase of accretion, during which they are observed as 
active galactic nuclei (AGN). AGN growth is intimately 
tied to galaxy evolution, as evident in the well-studied 
correlations between SMBH mas s {Mbh) and proper- 
ties of the h ost gala xy bii l ge fe.g.. |Geb hardt ct al.l |2000l : 
IFerrarese fc Merritti [2000t IMarconi "fc Hu nt 2003). The 
AGN phase is also hypothesized to regulate star forma- 
tion in its host g alaxy, with the galaxy feeding the black 
hole in turn (e.g.fPi Matteo. Springel. fc Hernquistil2005l : 
lYouneer et all 12008'). All massive galaxies are thought 
to experience ep isodic AGN beh avior in their lifetime 
(|Soltanlll98i iMarconi ei~alll200l . 

AGNs are generally classified by differences in their 
optical spectra. Type 1 or broad-line AGNs have broad 
{vpwHM ^ 1000 km s~^) emission lines superimposed 
on blue unobscured con tinua in the UV/optical (e.g., 
IVanden B erk et all [200ll ). and are the most luminous 
persistent sources in the sky. Type 2 or narrow-line 
AGNs lack broad emission lines and have weaker con- 
tinua (frequently dominated by their host galaxies), but 
have strong narrow emission lines, especially from for- 
bidden transitions. Narrow emission lines associated 
with nuclear activity can be distinguished from lines 
caused by star formation by study ing the line ratios 
(|Baldwin. Phihps. fc Terlevich|[T98l . The line ratio di- 
agnostics work because the "harder" emission of an 
AGN is more efficient at ionizing the surrounding gas 
and dust than star formation, and thus AGNs have 
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stronger lines from high-energy forbidden transitions 
(e.g., [Oiii] A5007A and [Nil] A6583A) relative to the 
lower-energy hydrogen transitions (e.g., H/3 A486lA and 
Ha A6563A). The subclass of "l ow-ionization n uclear 
emission region" AGNs (LINERs, lHeckman|[T980t ) have 
narrow emission lines that are probably excited by 
some combination of ionization frorn both star forma- 
tion and an AGN (Eracle ous et all 120101) . Deep X- 
ray sur veys have additionally revealed "optica lly dull" 
AGNs HElvis et all 119811: iComastri et al.l I2002D . which 
have bright X-ray emission but none of the broad 
or narrow emission line signatures of AGN accretion. 
While many optically dull AGNs can be explained 
as Type 2 AGNs diluted by prominent host galaxies 
(iMoran. Filippenko fc Chornockl [200^ iCaccianiga et al.l 
I2007D . at least ~l/3 are undiluted but intri nsically opti- 
cally weaker than other AGNs ([Trump et al.ll2009(J) . The 
inferred X-ray column density Nh can also be used to 
classify AGNs, with Type 2 (narrow-line) AGNs typically 
more X-ray absorbed than Type 1 (broad- line) AGNs. 
However X -ray and optical clas sifications differ for ^20% 
of objects (jTrouille et al.ll2009D . 

Historically, Type 2 and optically dull AGNs have 
been described as obscured versions of Type 1 AGN, 
with the broad emission line region (BLR) hidden be- 
hind a partially opaque "torus" of gas and dust, while 
the narrow emission line s lie outside the torus (e.g., 
iKrolik fc BegelmanI I1988D . The best evidence for this 
scenario is the observation that some Type 2 AGNs 
have a "hidden" BLR revealed by spectropolarimetry 
(jAntonuccil Il993| ). However, recent observations have 
revealed several serious limitations of a simple unified 
model based solely on geometric obscuration. Even in 
very deep spectropolarimetr ic observations, many Type 2 
AGN s show no hi d den BLR (iBartK Fihppenko fc MoranI 
IT999I : [TVMl[200lt I Wang fc Zhand 120071 ). Observations 
suggest a lower L/L Kdd > 0-01 limi t in accretio n rate 
for bro ad-line AGNs ([Kollmeier et al.li2006: Trump et al.l 
l2009bf) . although they rem ain incomplete at low accre- 
tion rates and low masses (IKellv et al.l[20T0l) . The X- 
ray spectra are unab sorbed (Nh 5, 10^^ cm~^) for 30- 
40% of Type 2 AGNs (jMainieri et a.l.|[2 007': 'TrouilleeLaD 
I2009D . as weh as most local LINERs (Ho 200|, and 
references therein) and distant optically dull AGNs 
(|Trump et al.ll2009c|) . Several well-studied LINERs ad- 
ditionally la ck the narrow Fe Kg emission signat ure of a 
dustv torus (|Ptak et al.ll2004l : lBianchi et al.ll2008l ). Many 
Type 2 AGNs and most opticall y dull AGNs have mid- 
IR co lors like normal galaxies ()Hol l2008t iTrump et al.l 
l2009c[ ). in contrast to hot mid-I R colors of Type 1 
AGNs described by torus mo dels (jNenkova et al.l l2008l : 
iMor. Netzer. fc ElitzuJ l2009| ) . Toroidal obscuration is 
a dditionally ruled out for some strongly varying Type 
2 ('Hawkinsl l260l and optically duU AGNs ([Trump et al.l 
|2009c), since these objects have continua which vary on 
year timescales, well within the inferred light travel time 
dimension of any torus. 

Several authors have proposed models which use differ- 
ent accretion rates a s a cause o f the differences between 
observed AGNs. Elit zur fc Hoi (2009) suggest that the 
BLR and "torus" are inner (ionized) and outer (clumpy 
and dusty) parts of the same disk-driven wind, and that 
this wind is no longer supported at low accretion rate 



(see also lEUtzur fc Shlosman| [2006l: iNenkova et "aIll2008D . 
Similarlv. iNicastroT i 2000D "suggested that low accretion 
rates actually drive the disk wind within the last sta- 
ble orbit of the SMBH, meaning that the BLR cannot 
form. Models for radiatively inefficient accretion (e.g., 
I YuanI [2OO7I ) suggest that at L/Lsdci £ 10~^, the accre- 
tion disk becomes truncated near the SMBH, with a geo- 
metrically thick and optic ally thin disk at low radii, and 
a normal thin disk (e.g., iShakura fc SunvaevI 119731 ) at 
higher radii. Such objects are predicted to lack strong 
emission lines (both broad and narrow) and have weak 
UV/optical emission, as observed in many optically weak 
low- luminosity AG Ns (IHoI 120091) and X - ray bright, op- 
tically duff AGNs (|Trump et al.ll2009c[ ). iHopkins et all 
([2009') additionally show that X-ray hardness, generally 
attributed to X-ray absorption, may also result from the 
naturally X-ray hard spectrum expected from radiatively 
inefficient accretion. 

In this work we directly measure Eddington ratios for 
a large. X-ray selected sample of broad-line, narrow- 
line, and lineless AGNs. The Eddington ratio is a 
unitless measure of accretion power, defined as A = 
Lint/ LEdd- (with Lint the intrinsic accretion luminos- 
ity). The sample is drawn from the Cosmic Evolution 
burvey fCOSMOS.IScoville et al.l l2007l) X -ray AGN sam- 
ple (jTrump et al.ll2007l ). as described in Section 2. Esti- 
mates of specific accretion rates are described in Section 
3, with intrinsic accretion luminosity Lint measured di- 
rectly from fits to the multiwavelength continuum (avoid- 
ing uncertain bolometric corrections) and black hole 
masses from the broad line scaling relations (for broad- 
line AGNs) or the Mbh — M^, relations (for narrow-line 
and lineless AGNs). In Section 4 we show that broad 
emission lines are present at only high accretion rates 
[Lint/ LEdd > 0.01), while narrow-line and lineless AGNs 
at lower accretion rates have cooler disks, stronger ra- 
dio jets, and no torus IR signature. We present a "car- 
toon" model which summarizes our results in Section 5, 
with predictions for future observations in Section 6. We 
adopt a cosmology with h = 0.70, JIm = 0.3, Q.a = 0.7 
throughout. 

2. OBSERVATIONAL DATA 

Measuring an accurate specific accretion rate re- 
quires accurate accretion luminosities and black hole 
mass estimates. In particular, SED measurements 
from optical/UV to X-ray are necessary to con- 
strain intrinsic luminosities to within a factor of 
a few (as we show in §3.1). We select a sam- 
ple of 348 AGNs from the Co smic Evolution Survey 
(COSMOS, IScoville et afl l2007l) field which is based 
on th e 1.7 deg^ HST/ACS mosaic (jKoekemoer et all 
|2007[ ). These AGNs have multiwavelength data in 
the form of Spitzer/IRAC, HST/ACS, Subaru/Suprime- 
Cam, GALEX, XMM-Newton, and Chandra observa- 
tions, as described in Table [TJ Spectroscopic identifica- 
tion and redshifts for these objects comes from archival 
SDSS data, Magellan/IMACS and MMT/Hectospec 
(Trump et al. 2009a), and VLT/VIMOS observations 
( Lilly et al.ll2007i) . 

The sample is selected from the parent catalog of 1651 
XMM-COSMOS p oint sources with optical counterparts 
(|Brusa et al.ll2010t ). limited by /o.5-2koV > 2 x IQ-^^ 
erg s~^ cm~^. Of these X-ray point sources, 649 objects 
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with ?AB < 23.5 have high-confidence (> 90% likelihood 
as correct ) identifications and redshifts from opt ical spec- 
troscopy ((Trump et al.ll2009at iLiUv et al.|[2007l ) in COS- 
MOS. Most of the X-ray point sources without spec- 
troscopy were missed simply due to random slit place- 
ment constraints. The optical spectroscopy is ^90% 
complete to iab < 22.5, although the completeness is 
redshift-dependent. For broad-line AGNs, the spectro- 
scopic completeness is lower at0.5<2<l, 
an d z ^ 2.4, e specially at zab > 22.5 (see Figure 13 
of iTrump etal . 2009aj. For narrow-line and lineless 
AGNs, spectroscopic completeness drops dramatically at 
z > 1.2, since at higher redshifts the 4000 A break and the 
[On] feature shift redward of the observed wavelength 
range. To ensure that X-ray objects with narrow-line and 
lineless spectra are bona-fide AGNs, we select only ob- 
jects with Lo.5-iokcV > 3x 10**^ erg s~^. This X-ray lumi- 
nosity limit is generally used to se parate AGNs from X- 
ray f ainter starburst galaxies fe.g.. iHornschemeier et ahl 
[200l . We also include seven broad-line AGNs with- 
out X-ray detection, six of which were selected by their 
Spitzer/IRAC colors and one which is a serendipitous ob- 
ject from the bright zCOSMOS survey (which selected 
targets based only on iab < 22.5). While these 7 X-ray 
undetected AGNs do not come from a complete sample, 
we include them to gain a larger parameter space of AGN 
spectral types and accretion rates (in effect, when using 
their X-ray limits, they occupy the same Ldisk/Lx pa- 
rameter space as a few other X-ray detected AGNs). Re- 
stricting narrow-line and lineless AGNs to be X-ray lumi- 
nous and adding the 7 X-ray undetected broad-line AGNs 
makes a parent sample of 380 broad-line, 124 narrow-line, 
and 49 lineless AGNs (553 total) with high-confidence 
redshifts and spectral identification. 

Measuring accurate black hole masses additionally con- 
strains the sample to certain redshift ranges. For Type 1 
AGNs, we require the presence of one of the Civ, Mgii, 
or H/3 broad emission lines in the observed spectral range, 
effectively limiting broad-line AGNs with IMACS or VI- 
MOS spectra to 0.16 < z < 0.88, I < z < 2.4, and 
2.7 < 2; < 4.9 and objects with Hectospec or SDSS spec- 
tra to z < 4.9. For narrow-line and lineless AGNs, we es- 
timate black hole mass from the Mbh ~ L^uige relation, 
and so we require an accurate estimate of Lbuige- For 
this we use the sample of o bjects in COSMOS with mor- 
phological decompositi ons (|Gabor et al.1 12 0091 from the 
HST/ACS images (Koekemoer et al. 020071 which also 
effectively limits the narrow-line and lineless AGNs to 
z < 1.2 (beyond which the 4000 A break shifts out of 
the ACS-i band and the host galaxy is much more dif- 
ficult to detect). T he accurate host measurements from 
iGabor et al.l (|2009[ ) additionally allow us to subtract the 
host component before computing the intrinsic bolomet- 
ric luminosity. In general, the narrow-line and lineless 
AGNs are biased towards lower redshift and consequently 
higher mass, since AGNs grow over the cosmic time. The 
narrow-line and lineless AGNs have a mean redshift of 
0.7, while the broad-line AGNs have a mean redshift of 
1.6. The final sample of 348 AGNs includes 256 broad- 
line, 65 narrow-line, and 27 lineless AGNs. 

Full multiwavelength data exist for > 95% of the AGNs 
in the sample in every wavelength region except the 
UV. X-ray data exist from both Chandra and XMM- 



Newton: we use the deeper Chandra data when available, 
but the Chandra observations cover only the central 0.8 
deg2 of the COSMOS field. For the 7 X-ray undetected 
broad-line AGNs, we use the 0.5-2 keV XMM flux limit 
(/o.5-2kcV — 2 X 10~^^ erg s~^ cm~^) for their X-ray lu- 
minosity (since these AGNs have Ldisk/Lx > 10, their 
bolometric luminosity is dominated by their optical/UV 
emission and completely neglecting their X-ray emission 
does not significantly change their bolometric luminos- 
ity estimate ) . W e apply the zero-point offsets derived by 
lllbert et al.l doM) to the IR-UV photometry. 

2.1. Measuring Absorption and Extinction 

X-ray absorption and optical/UV extinction could pose 
a challenge to measuring the intrinsic accretion power. 
The most heavily absorbed AGNs (e.g. Compton-thick 
AGNs with Nh > 10^^ cm~^) are entirely missed by our 
survey because they lack detectable X-ray emission (e.g., 
iTreisteret "all 120041 ). But if an AGN is moderately ab- 
sorbed and still X-ray detected, we might expect its disk 
to appear cooler because the UV light is preferentially ex- 
tincted, and its X-ray slope to appear harder because the 
soft X-rays are preferentially absorbed. Some AGNs are 
also intrinsically reddened, decreasing their UV emission 
by a factor of 2-3 (Richards ct al. 2003) and causing us 
to underestimate their accretion disk emission. With ab- 
sorbed soft X-rays and extincted disk emission, we could 
significantly underestimate Lint/ LEdd- 

We use X-ray column density Nh to characterize 
the obscuration properties of our AGNs. Column den- 
sity and optical extinction arc rough ly correlated, wit h 
Ay /Nh --2x10^23 cm^ (Martinez-Sansigre et al.ll2006D . 
Then at Nh < 10^^ cm~^, optical magnitude should 
be extincted by < 20% (< 0.2 mag). Assuming a 
SMC reddening law (|Peilll992D . as is most appropriate 
for AGNs, this optical extinction translates to a factor 
of ^1 2 extinction at 3000Ain the UV. iMaiohno erahl 
(|2001f ) showed that the Ay — Nh relation varies by up to 
a factor of 30 because of unknown changes in the gas-to- 
dust ratio, grain size, and/or different physical locations 
of the optical and X-r ay absorbing material. However 
for all AGNs in the .Maiolino et al.l ([2001) sample with 
Lx > iC^ erg s~^, Ay/Nn < 1.8 x 10"^^ cm^, meaning 
at Nh ~ 10^^ cm~^ even the maximum optical (y-band) 
extinction is a factor of 5 and the maximum UV (3000A) 
extinction a factor of 30. 

Column density Nh can be accurately measured for 
the 153 AGNs (93 broad-line, 38 narrow-line, 22 line- 
less AGNs) in the sample which have > 40 XMM or 
Chandra counts. (With less than 40 counts, the spec- 
tral fitting does not always stably converge.) We fit each 
X-ray spectrum as an intrinsically absorbed power-law 
with Galactic absorption {Nn^gai = 2.6 x 10^° cm^ in 
the direction of the COSMOS field), with the power- law 
slope and Nh as free parameters. The best-fit Nh value 
and its error are found using the ICashl (fT979h statistic. 
We present Nh and X-ray slope F in Figure [TJ Among 
the 153 AGNs with > 40 X-ray counts, there are 118 
unobscured AGNs with Nh < 10^^ cm"^ (82 broad-line, 
24 narrow-line, and 12 lineless AGNs). We restrict our 
main conclusions to this set of 118 unobscured AGNs for 
the remainder of this work. 
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Fig. 1. — The column density Nh and the X-ray slope Tx mea- 
sured from the X-ray spectrum for the 153 AGNs with > 40 XMM 
or Chandra counts. X-ray slope Fx is defined by L 
Black crosses show broad-line AGNs, blue diamonds show narrow- 
line AGNs, and red squares show lineless (optically dull) AGNs. 
The median X-ray slope for all AGNs is Vx = 2.1, although 
Vx ranges from 1 to 3. There are 118 unobscured AGNs with 
Nh < 10^2 cm-2. 

3. CHARACTERIZING AGN SPECIFIC ACCRETION RATE 

In this work we describe the specific accretion rate us- 
ing the Eddington ratio parameter, A = Lint/ LEdd- Here 
Lint is the intrinsic luminosity, a measure of the total 
accretion luminosity which includes only light from the 
accretion disk and X-ray corona and excludes any repro- 
cessed IR emission. While the reprocessed IR emission 
can represent a large fraction of the bolometric luminos- 
ity, especially for obscured AGNs, it may be anisotropic. 
Most of our AGNs are unobscured (see Section 2.1) and 
we exclude the IR emission to avoid double-counting 



the AGN emission. Instead we use only the optical/UV 
and X-ray emission which comes directly from the disk 
and corona in the AGN: in this work when using "in- 
trinsic" luminosity we are always referring to the total 
of the disk (optical/UV) and corona (X-ray) emission, 
without the reprocessed (IR) emission. The Eddington 
luminosity is derived from the black hole mass, with 
LEdd = 1-3 X W^{Mbh/Mq) erg 8"^ AGN luminos- 
ity is powered by accretion rate, with Lint = rjMc^. 
For a constant efficiency rj, the Eddington ratio A is 
equivalent to the specific accretion rate m = M/MEdd- 
For example, assuming 77 0.1 the Eddington accre- 
tion rate can be written MEdd = 5Mg Mq yr~^ with 
Ms = M/{10^Mq). However there is good evidence 
that 77 decreases at very low accre tion rates m << 0.01 
(e.g.. iNaravan fc McClintockl [20081 ) . Indeed, in Sections 
4 and 5 we invoke a lower- efficiency (radiatively ineffi- 
cient) accretion to explain the observational properties 
of Lint/ LEdd < 10~^ AGNs. This means that the accre- 
tion power Lint/ LEdd probably underestimates the ac- 
cretion rate m for our most weakly accreting AGNs with 
Lint/ LEdd < 10~^: e.g., a measured accretion power of 
Lint/LEdd ~ 10"'' might correspond to m ^ 10"'^. 

Below we outline our methods for estimating black hole 
masses and bolometric luminosities from the data for the 
AGNs in our sample. Table [2] presents the full catalog of 
Lint, Mbh, and Lint/LEdd, and their associated errors, 
for our AGNs. 

3.1. Intrinsic Luminosity Estimates 

We calculate the intrinsic luminosity from the full rest- 
frame near-IR to X-ray multiwavelength data. This 
avoids monochromatic bolometric corrections which are 
highly un certain and probably depend on Eddington ra- 
tio fe.g., iKellv eral...2008; .Vasudevan fc Fabian .2009.) . 
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Instead we measure intrinsic luminosity by integrating 
the best-fit accretion disk -I- X-ray power-law SED model. 
We compile the broad-band near-IR {Ks, J), optical {z'^ , 
r+, i\ 5+, Vj, Bj, u*), UV (GALEX NUV & FUV), and 
X-ray (0.5-2 keV and 2-10 keV from Chandra when avail- 
able or XMM-Newton) data, for which the wavebands 
and limits are described in Table [T] To avoid repro- 
cessed mid-IR emission, which would double-count the 
intrinsic emission for an unobscured AGN, we restrict 
the accretion disk fit to rest-frame 1 < i? < 100 eV 
(6200 > A > 124 A). The radio emission is negligible 
in the total energy output of our AGNs. While narrow- 
band optical photometry also exists for our AGNs, its 
inclusion doesn't appreciably change the best-fit multi- 
wavelength SED from using only the broad-band data. 

The rest-frame near-IR and optical emission of narrow- 
line and lineless AGNs is dominated by the emission 
from the host galaxy. For these objects, accurate in- 
trinsic luminositie s require modeling and subtracting the 
host galaxy light. iGabor et al.l (|2009f ) measured the host 
F814W luminosities from surface brightness fitting to 
the HST/ACS data of our AGNs. We use thi s lumi- 
nosity to scale a galaxy template from PoUet ta et al.l 
(j2007l) . Lineless AGNs have early- type hosts, since 
their spectra lack the emission lines associated with 
a late- type star- forming galaxy, and so we u se th e 
"E115" early- type template from IPolletta et all (|2007[ ). 
The narrow-line galaxies in our sample typically have 
intermediate- t ype ("green v alley") hosts based on their 
morph ologies (|Gabor et al? 2009) and star formation 



rates (iSilverman et al.l 1 2009). and so we use the "SO" 
template of IPolletta et al., (|2007i) . We subtract the host 
contribution in each photometric band before perform- 
ing our SED fit. The redde st ("E112") and blues t ("Sd") 
normal galaxy templates of IPolletta et al.l (120071 ) are ad- 
ditionally used as extreme hosts to estimate the possible 
error contribution from choosing the wrong host tem- 
plate (described in §3.3). 

It is possible that a few of the narrow-line and line- 
less AGNs might have very blue starbursting hosts, al- 
though such galaxies are uncommon at z < 1. An ex- 
tremely blue, UV-emitting host would cause us to over- 
estimate the accretion disk emission and consequently 
overestimate the accretion rate. Since the narrow-line 
and lineless AGNs have lower accretion rates than broad- 
line AGNs (as we discuss in Section 4), if their true ac- 
cretion rates were even lower it would only strengthen 
our conclusions. It is also possible that very red, dusty 
hosts could cause us to underestimate the true accretion 
rates for narrow-line and lineless AGNs. However a dusty 
host should cause the AGN to appear extincted, and our 
sample of AGNs generally has low measured absorption 
(see Section 2.1). In addition, restricting our fitting to 
1 < < 100 eV (6200 > A > 124 A) aheady means that 
a normal elliptical galaxy (like our "E112" template) con- 
tributes very little flux where we fit the accretion disk. 

While broad-line AGNs are likely to have some host 
contribution, we cannot use surface brightness fitting 
to estimate their host luminosity because they are at 
high redshift and their point source o verwhelms their ex- 
tended emission (jGabor et al.l I2OO90 . However at the 
peak of the accretion disk emission for a broad-line 
AGN (-- 3000A, or 4 eV) the host galaxy contributes 



< 20% of the emission fe.g.. iBentz etalll2006l ). Be- 
cause we additionally restrict our accretion disk fitting 
to 1 < < 100 eV (6200 > A > 124 A), we can as- 
sume that the error from not subtracting the host for 
broad-line AGNs is typically < O.ldex. 

We shift the observed (and host-subtracted, for 
narrow-line and lineless AGNs) photometry to the rest- 
frame from the measured spectroscopic redshift and con- 
vert the magnitudes or fluxes to luminosities. We then 
fit an accretion disk model to the optical/UV emission 
within the range 1 < < 100 eV (4.8 x lO" < p < 
2.4 X 10^^ Hz, or 6200 > A > 124 A) and a power-law 
representing the X-ray corona emission to the rest-frame 
X-ray data. We measure the total bolometric luminosity 
from the sum of the disk luminosity (given by the an- 
alytic solution in Equation 3 below) and the power-law 
luminosity from AEpeak < E < 250 keV (where Ep^ak is 
the peak energy of the best-fit disk model) . While the X- 
ray background requires a h igh-energy cutoff for AGNs 
in the few hundreds of keV (jGilli et al.ll2007[) . measure- 
ments of the cutoff ene rgy exist for only ^ 15 AGNs and 
vary from 50-500 keV (jPerola et al.l[200l iMolina et all 
I200I . We choose 250 keV as an intermediate value, 
although any cutoff from 50-500 keV does not greatly 
influence our results. Our AGNs have typically flat X- 
ray spectra with Tx 2, and so changing the X-ray 
cutoff energy by a factor of 0.2-2 effectively changes the 
integrated X-ray luminosity by the same factor of a few. 
Because the X-ray and disk luminosities are roughly com- 
parable (see Figure [2]), this results in less than a factor of 
two change in the total accretion luminosity: much less 
than the '^0.5 dex errors we compute for our estimated 
Lmt/LEdd (see §3.3). ^ 

We use the accretion disk model of IGierlihski et al.l 
Ip99l ). which improves upon a basic blackbody accre- 
tion disk by including a cor rection for relativistic effects. 
fThe lGierlihski et all (|1999^ model is the "diskpn" model 
of the XSpec X-ray fitting software.) This model is 
based on the pseudo-Newtonian gra vitational p otential 
$ = -GM/{R - Rg) (iPacvnski fc Wlitai [T98ol) . where 
R„ is the Sch\yarzsch ild radius Rg = GM/c^. From 
IGierlihski et al.l ()1999f ). the model takes the form: 



L = KE'' 



rdr 



exp[E/kT{r)] - 1 



(1) 



where r — R/ Rg and we assume the innermost stable 
orbit Tin — 6. The temperature depends on radius as 



Co 



2/3 



r(r - 2)3 



1 - 



33/2(^ 



2l/2y,3/2 



1/4 



(2) 



with Co ~ 0.1067, and To oc m^^/^m^/*. The coefficient 
K depends on inclination angle, coronal absorption, and 
the color to effective temperature ratio. Rather than 
estimate these values, we assume that X is a constant, 
computed by simply scaling the model to our data. To is 
the sole free parameter. In our analyses below we refer to 
Epeak, the peak energy of the disk, rather than Tq, and in 
general kTo ~ Epeak/'^^:. We find the best-fit disk model 
in terms of Tq by minimizing the function. While 
most of the best-fit disk models have significant emission 
at E < 1 eV, we restrict the fit to 1 < i? < 100 eV to 
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mitigate the effects of a contaminating torus and/or host 
galaxy light. 

Note that the relation Tq oc m~^/'*m^/'' above means 
that the disk temperature is constrained not only by the 
photometry by also by the black hole mass. In prac- 
tice this prevents our fits from resulting in unphysically 
hot accretion disks, since disks peaking at energies much 
higher than -4 eV (3000A) would require unphysically 
small black hole masses. This is especially important to 
note because about one-third of the sample lacks G ALEX 
UV detections, and as a result the declining high-energy 
slope is not well constrained by the photometry for low 
redshift AGNs. The black hole mass error (~ 0.4 dex) is 
used during the bootstrapped uncertainty measurements 
for the accretion disk temperature and luminosity. 

The total dis k luminosity is ca.lculat ed analytically (see 
Appendix A of iGierlihski et al.l[T999( ): 



-'disk 



K 



167r 



Co 



(3) 



Errors in both Ep^ak and Ldisk are found by bootstrap- 
ping 1000 fits to the resampled data. 

To characterize the X-ray corona emission, we use the 
X-ray spectral fits described in §2.1. Each X-ray spec- 
trum is fit as an intrinsically absorbed power-law with 
Galactic absorption {Nn.ged — 2.6 x 10^" cm^ in the direc- 
tion of the COSMOS field). We use the photon index Fx 
to represent the power-law slope, such that L^, = Lqv^~^ . 
Figure [1] shows that the typical Tx — 1.9±0.4, and we as- 
sume this slope for AGNs with too few X-ray counts for a 
good fit. We calculate the total X-ray luminosity by inte- 
grating the power-law model over AEp^ak < E < 250 keV 
(where Epeak is the energy peak of the disk model) , using 
the analytic solution: 



Lx = io/(2-r)x[(250keV//i)^-i -{iEpeak/h) 



\2-ri 



(4) 



The total bolometric luminosity is simply the sum of 
the integrated accretion disk and X-ray power-law com- 
ponents. Lint = Ldisk + Lx- 

Figure [5] shows a representative sample of broad-line, 
narrow-line, and lineless SEDs with model fits. Note 
that emission lines and variability (the various photo- 
metric data were taken over 3 years) mean that our sim- 
ple accretion disk model is not a perfect fit: some of the 
optical/UV data and differ from the model fit by up to 
0.2 dex. However, such small errors in individual pho- 
tometry points are negligible compared to the > 0.4 dex 
total errors we estimate for Lint (see Figure|4]and Section 
3.3). In general, the accretion disk plus X-ray power-law 
model provides an accurate, physically motivated fit to 
the data. 

3.2. Black Hole Mass Estimates 

For Type 1 AGNs, we estim ate black hole masses 
using the scaling relations of IVestergaard fc Osmeri 
(I2009f) for the Mgi i broad emission line and 
IVestergaard fc PetersonI (|2006n for the H/3 and 
Civ broad emission lines. These relations estimate 
black hole mass from single-epoch spectra by employ- 
ing the correlation between the radius of the broad 
emission line region and the continuum luminosity. 



Rblr ^ L '^'^, observed in loc a l AGN with r everb eration 
mapping ijBentz et al.l 120061 : iKaspi et al.l 120071 ). In 
general, masses estimate d from the scaling relations are 
accurate to ^ 0.4 dex (jVestergaard &: PetersonI 120061 : 
IShen e t aL 2008) and agree with local AGN masses from 
dyna mical estimators (Da vies et al.l 20061: lOnken et ahl 
20071) and the Mt^ h-u* correlation (jOnken et al.l 12004 



Greene fc IIall2066l ). The scaling relations take the form 



of Equation 5, with \Lx in units of 10** erg/s and 
vpwHM in units of 1000 km/s; A = 6.91, B — 0.50, 
and A = 5100A for H/3; A = 6.86, B = 0.50, and 
A = 3000A for Mgii; A = 6.66, B = 0.53, and 
A = 1350 A for Civ. 



log 



BH 



V Mr, 



= A + B \og{XLx) + 2 \og{vFWHM) (5) 



Black hole masses for the Type 1 AGNs with Magel- 
lan/IMACS or SDSS spectra in COSMOS are a lready 
published in previous work (| Trump et al.l l2009b( ) , and 
we repeat the same techniques for Type 1 AGNs with 
VLT/VIMOS spectra. Briefly, a power-law fit plus iron 
emission are fit to each AGN. The continuum luminosity 
is estimated directly from the continuum fit, while the 
velocity widths are computed from Gaussian fits to the 
continuum-subtracted emissio n lines. Some o bjects also 
have black hole masses from iMerloni et al.l (|2010( ) ; for 
these objects, our masses are consistent with a random 
scatter of only ~0.4 dex: equivalent to the intrinsic scat- 
ter of the scaling relatio n s (see Figure 3 of iTrump et al.l 
I2009bf l. iMarconi et all (|2008D showed that the scat- 
ter in Mbh from the scaling relations might decrease 
to 0.2 dex if radiation pressure is taken into account. 
Replacing the scaling relations from Equation 5 with 
those of ,Marconi et al.l p008. ) would tighten the distribu- 
tion of Lint/LEdd estimates for broad-line AGNs about 
Lint/L^dd ^ 0.3. This has no impact on the Mbh es- 
timates for narrow-line and lineless AGNs, and does not 
affect the difference in Lint/ L^dd between the broad- line 
sample and the narrow-line and lineless AGN sample. 

Estimating black hole masses for AGNs without broad 
emission lines requires secondary estimators. We employ 
the relationship betwe en Mb h and rest-frame X-band 
host bulge luminosity (jGrahamll2007| ): 



log 



M 



BH 



Ma 



0.93(log(iK)-0.3z)- 32.30, (6) 



with Lk in units of erg s~^. The Mbh — Lk, bulge relation 
comes from the more fundamental Mbh — Af* relation, 
since res t-frame K bulge luminosity is correlated with 
M* (e.g.. lllbert et al.ll2010l) . We add an additional -0.3z 
term to the ( Grahan] l2007f l relation in order to account 
for the evolution in the M^ /Lk ratio, \og{M^/ Lk) oc 
—0.3^ (|Arnouts et al.l I2007D . We measure rest-frame 
Lk from the host galaxy template from the multiwave- 
length SED fit (described above in §3.1). The early-type 
template for the lineless AGNs is, by definition, bulge- 
dominated, and so Lk. bulge = LK,host- The SO template 
used for the narrow-line AGNs, however, has a significant 
disk component, and so we take LK,buige = Q-^LK.host- 
The i ntrinsic error in the Mbh — Lk is 0.35 dex (Grahaml 
'20071). We do not correct the Mbh estimates for any 
evolution in the Mbuige — Mbh relation because mea- 
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Fig. 2. — Multiwavelength photometry and model fits for 12 example AGNs. The top four panels are broad-line AGNs (represented by 
'BL'), the middle four are narrow-line AGNs (represented by 'NL'), and the bottom four are lineless "optically dull" AGNs (represented 
by 'OD'). In each panel, the dashed line is the best-fit accretion disk model and the dot-dashed line is the X-ray power-law fit. The X-ray 
power-law slope comes from the X-ray spectral fit, although we show only the X-ray photometry data in this figure. Estimated host SEDs 
are shown by solid lines for the narrow-line and lineless AGNs. We fit only at > 1 keV in order to ignore the reprocessed IR emission, 
and so the longest- wavelength photometry data (especially the IRAC channels) are not fit by our models. 



suring of Mbuige — Mbh evolution h as proved difficult 
due to significant bi ases in most tests (jLauer et al.ll2007t 
IShen fc Kellvl l201^. Besides while ther e IS some evi- 
dence for evolution to ^ ~ 3 (jPecarli et a l. 2010 ), there 
is pro bably little or no evolution to z 1.5 (|Jahnke et alj 
mM) and' our narrow-line and lineless AGNs lie at < 1 . 

Because we use different mass estimators for broad-line 
and narrow- line/lineless AGNs, it is important to demon- 
strate that the two methods agree. Seven of our broad- 
line AGNs have detected host galaxies from the decom- 
positions of iGabor et al.l (|2009f) . and for these AGNs we 
compare Mbh estimates from the broad-line scaling re- 
lations and from the host galaxy rest-frame Lk in Figure 
[21 The Mbh estimates from broad lines and agree 
within < 2tT for all objects (indeed, estimates for all ob- 
jects but one agree within < Icr). In addition to the seven 
broad-line AGNs in our sample, both the broad-line and 
host galaxy Mbh estimators have bee n shown to produce 
consistent masses for nearby AGNs (jOnken et al.ll206^ 
iGreene fc Hol[2006[ ) . It is particularly unlikely that either 
of the estimators is systematically off by a factor of 100. 



Therefore we are confident that the factor of 100 differ- 
ence in Lint/ LEdd for broad- line and narrow-line/lineless 
AGNs in Section 4 (see, for example. Figure [5]) is a phys- 
ical effect, robust beyond the choice of black hole mass 
estimator. 

We highlight the range and limitations of the AGN 
sample in Figure El which shows bolometric luminosities 
and black hole masses for the broad-line, narrow-line, 
and lineless AGNs. Objects in the upper left have the 
highest specific accretion rates, while those in the lower 
right are weakly accreting AGNs. While the total sam- 
ple spans 3 orders of magnitude in both luminosity and 
black hole mass, our narrow-line and lineless AGNs are 
generally less luminous and more massive than broad- 
line AGNs. The lack of low-mass narrow-line and line- 
less AGNs is due to the selection limits of the survey: 
such objects are too faint to be detected in COSMOS. 
It is suggestive that these higher mass narrow-line and 
lineless AGNs are at z < 1 and are less luminous: this is 
consistent with " downsizing," with m ore massive AGNs 
becoming less active at lower redshift ()Ueda et al.ll2003l : 
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Fig. 3. — Black hole mass estimates from both the host Lk and 
the broad-line scaling relations for the sev en broad-line AGNs with 
detected host galaxies from G abor et al.l ((2)09). For all but one 
AGN, both Mgjj estimates agree within < Icr (for the remaining 
object, the two estimates differ by only ~ 2a). From these AGNs, 
and the sets of nearby AGNs with similarly consistent masses from 
both estimators (Onken et al. 2004; Greene & Ho 2003), it is un- 
likely that the different mass estimators cause bias between broad- 
line and narrow-line/lineless AGNs. 
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Fig. 4. — Intrinsic luminosity Li„t with black hole mass Mbh 
for the AGN sample. Broad-line AGNs are shown by black crosses, 
narrow-line AGNs by blue diamonds, and lineless AGNs by red 
squares. Errors are calculated as described in Section 3.3. Narrow- 
line and lineless AGNs generally have higher mass, due to the 
COSMOS selection limits, but they also have lower luminosities 
as expected by downsizing. 

iBrandt fc Hasingeii[2005l: iBongiorno et al.ll2007[ ). 

Figure 3] shows that at a given mass or luminosity there 
are generally all types of AGNs present in our sample. 
For this reason we do not expect that the differences 
between broad- line and narrow-line/lineless AGNs are 
biased by selected samples from different masses or lu- 
minosities. In addition, despite the different redshifts 
of most broad-line and narrow-line/lineless AGNs, we 
do not expect their differences to be caused by red- 
shift. There is evidenc e that AGN obscuration proper- 
ties de pend on redshift ([Treister et al.|[2009l : lTrump et al.l 
I2009al) . but these AGNs are unobscured. The AGN cen- 
tral engine, meanwhile, does not change with redshift 
in terms of ionization parameters ([Dietrich Sz HamannI 



I2004t lVestergaardll2004f). spectral e nergy distributions 
(Vignah et al.l 120031 iRichards et al.l [2001 iKellv et al.l 
12008) . or metahicity ()Simon et al.l I2O1O0 . Limiting the 
sample to 2: < 1, 8.5 < log(MBjy) < 9, or 44 < 
log(Lmt) < 45 does not significantly change the dif- 
ferences between the broad-line and narrow-line/lineless 
AGN samples seen in Figures 5, 6, 7, or 8. 

3.3. Error Budget 

We estimate errors for each of our specific accretion 
rates, propagating the errors from both the intrinsic lu- 
minosity estimate and the black hole mass estimate. Our 
intrinsic luminosity is subject to three major uncertain- 
ties: 

• Photometry errors, Uphot- We measure the er- 
ror contribution of the photometry by bootstrap- 
ping, fitting our model SED to 1000 realizations 
of randomly drawn photometry values distributed 
according to the measurement errors. In general, 
o-phot ~ 0.1 dex. 

• Errors in the host subtraction, ahost- For broad- 
line AGN, we do not subtract a host component 
and assume that any remaining galaxy light overes- 
timates the intrinsic luminosity (from the UV and 
X-ray) by only < O.ldex (see §3.1). For narrow- 
line and lineless AGN we estimate ahost from the 
difference in the resultant Lint when using a very 
red ("E112 " ) and a very blue ( "Sd" ) template from 
lPolletta~et all (|2007i ). Since the accretion disk is fit 
only a.t E > 1 eV where there is little host emission 
(even from the "Sd" galaxy), this error is usually 
insignificant {ahost ^0.1 dex). 

• Incorrect Lint resulting from extinction, aext- Ex- 
tinction will make the true Ldisk greater than our 
estimate because optical/UV light will be missed, 
but will make the true Lx lower than our estimate 
because the power-law slope will be too hard. Be- 
cause we restrict our main analyses to unobscured 
{Nh < 10^^ cm~^) AGNs, we assume this error is 
< 0.1 dex (see §2.1). 

The black hole estimate is subject to two major uncer- 
tainties: 

• Intrinsic errors in the Mbh relations, arei- For 
broad-line AGN, the int rinsic error in the scal- 
ing r elations is 0.4 dex ([Vestergaard fc PetersonI 
I2006D . such that arei = 2.5Mbh- For narrow-line 
and lineless AGN, we use the Mbh — Lx.host rela- 
tion, and its associated intrinsic scatter is 0.35 dex 
(jGrahamll200'7t l. such that arei/MBH = 2.2. These 
errors dominate the error in Lint/ LecIcI, except for 
highly absorbed AGNs with Nh > 10^^-^ cm^^. 

• Measurement error in the luminosity used in the 
scaling relation, aium- For broad- line AGN, this 
is the measured continuum luminosity associated 
with the appr opriate scaling relation, estimated by 
iTrump et al.l (|2009b,) as aium 0.05 dex. Since 
Mbh oc L"-^, aium = I.SMbh for broad-line 
AGNs. For other AGNs the aium comes from our 
measured Lx^rest- We estimate this error for the 
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Fig. 5. — The distribution of calculated specific accretion rates 
(Lint/LEdd), for the 82 unobscured (A^j^ < lO'^^ cm~^) broad- 
line AGNs (black histogram), 24 narrow-line AGNs (blue dashed 
histogram), and 12 lineless AGNs (red dotted histogram). Narrow- 
line and lineless AGNs have significantly lower accretion rates than 
broad-line AGNs. The Lint/ L^dd 0-01 limit for broad-line 
AGNs is not a selection efl'ect IITrump et al.ll2009bl ). 

narrow-line and lineless AGNs from 1000 fits to the 
randomly subsampled data, and find that the er- 
ror is generally insignificant compared to the in- 
trinsic error {aium ^ 0.05 dex). Note the contri- 
bution from error in vpwHM to Mbh in broad- 
line AGNs is also negligi ble, since for our AG Ns 
(t{vfwhm) < 0.2VFWHM (iTrump et a ll l2009b( ). 

The total error in specific accretion rate, cta, is then 
given by: 



A2 



phot 



host 



7^2 
^int 



(7) 



We measure the total error a\ by bootstrapping, with 
1000 fits to the resampled data. In each fit we allow 
all of the parameters above to vary according to their 
error. The intrinsic error in the Mbh relations {urei) 
dominates the error. The average errors are ~0.5 dex, 
compared to the '^4 dex range in Lint / L Edd for the AGN 
in the sample. 

4. THE PHYSICAL EFFECTS OF SPECIFIC ACCRETION 

RATE 

The distribution of Lint/ L Edd for the 118 unobscured 
AGNs is shown in Figure [SI It is immediately evident 
that unobscured narrow-line and lineless AGNs accrete 
much more weakly than broad-line AGNs, with specific 
accretion rates differing, on average, by ^2 orders of 
magnitude. This suggests that many narrow-line and 
lineless AGNs are not simply geometrically obscured ver- 
sions of broad-line AGNs, but they instead have funda- 
mentally different accretion physics which we examine in 
more detail below. 

The large ^^0.5 dex errors in accretion rate artificially 
broaden the distributions, such that the intrinsic distri- 
butions are likely to be narrower than the histograms 
in Figure [5] appear (although many Lint/Lsdd ^ 1Q~^ 
narrow-line and lineless AGNs could be too faint for 
the COSMOS X-ray and spectroscopy limits). The 



Lint/LEdd ^ 0.01 limit for broad-line AGNs could 
be p artially explained by selection effects (|Kellv et al.l 
since low accretion rates AGNs are typically less 
luminous. However at the highest masses {Mbh ~ 
lO^M©), broad-line AG Ns with Lj^t/LEdd 5, 0.01 mus t 
be very rare (KoUmeier et al.ll2006l: [Trump et al.|[2009bl) . 
Meanwhile unobscured narrow-line and lineless AGNs 
are generally limited by Lint/LEdd ^ 0.01. With low 
X-ray column densities and low accretion rates, these 
objects hav e simi lar properties to the "naked" Type 
2 AGNs of iTranI ((2003), which additionally lack re- 
flect ed_broad_emission lines^^ (see 
also IGliozzi eFaII[2007t IWang fc Zhanl [20071) . We ex- 
pect that the X-ray unobscured low accretion rate AGNs 
would similarly lack reflected broad emission lines. Our 
method cannot accurately estimate Lint/LEdd for ob- 
scured AGNs, but fol lowing a un ifled model with geomet- 
ric obscuration (e.g.. lAntonuccil [7993.) . obscured narrow- 
line AGNs would likely have accretion rates comparable 
to our broad- line AGNs. 

We can compare the speciflc accretion rates and AGN 
types with the physical parameters of our model fits: 
namely the ratio of disk to power-law emission, the 
peak energy of the accretion disk model, and the X- 
ray power-law slope. These quantities are particu- 
larly useful in unifying AGN in terms of their accre- 
tion physics. Figure [H shows the specific accretion rate 
with these parameters for each AGN type. The val- 
ues of Ldisk/Lx can be roughly translated to values 
of gpy, with anx = - 3841og[L^(2500A )/L„(2keV)1 
([Tananbaum et al.l [19791 : iKellv et alJ [20081 ). The left 
panel of Figure [6| shows tracks of aox = 1.0,1.5,2.0, 
assuming Ep^ak — 6 keV and Tx = 1.9 (hotter disks and 
softer X-ray slopes increase aox)- Once again, narrow- 
line and lineless AGNs have lower specific accretion rates, 
and they also tend to have lower Ldisk / Lx and Ep^ak ■ 

We can determine the significance of any differences 
in Ldisk/Lx, Epeak, and Tx between rapidly accreting 
broad-line AGNs and weakly accreting narrow-line and 
lineless AGNs by comparing their mean values and con- 
sidering the scatter of each sample. Given mean val- 
ues iJi and /i2 and associated scatters Ci and (T2 for 
each set, the significance of their difference is given by 
(Ail - ^^2)/^f{(rl/Nl + crI/iVa) (where A^i and N2 are 
the numbers of AGNs in each sample. The broad-line 
AGNs have n(\og{Ld^sk/Lx)) = -0.14 ± 0.44 while the 
narrow-line and lineless AGNs have ^{\og{Ldisk/ Lx)) = 
—0.38 ±0.64, so that their difference is marginally signifi- 
cant at 2.1(7. The difference in Ep^ak is more significant: 
the broad-line AGNs have ti(\og{Epeak)) = 0.80 ± 0.20 
and the narrow-line/lineless AGNs have n{\og{Epeak)) = 
0.59 ± 0.37, so that the difference is significant to 3.3cr. 
From this we can conclude that a transition from weakly 
accreting narrow-line and lineless AGNs to the rapidly 
accreting broad-line AGNs results in significantly hotter 
and marginally brighter emission from the accretion disk. 

There is no significant difference between X-ray slope 
Tx for the different AGN types: mean Tx = 2.14 ± 0.29 
for rapidly accreting broad-line AGNs, and mean Tx — 
2.05 ± 0.29 for weakly accreting narrow-line and lineless 
AGNs (the difference is only 1.2(t sig nificant ). This is in 
contrast to the prediction of Hopki ns et al.l (|2009i) . who 
suggest that harder X-ray slopes are expected for radia- 




Fig. 6. — Specific accretion rate Li„t / LEdd ^.nd the ratio of disk to corona emission log{LiHgf;/Lx), disk temperature Epg^k and X-ray 
photon index Fx for the 118 unobscured AGNs with Nu < lO'^^ cm~'^. In each panel, black crosses represent broad-line AGNs, blue 
diamonds are narrow-line AGNs, and red squares are lineless AGNs. The dashed lines in the left panel show lines of apx = 1-0, 1.5, 2.0, 
assuming Ep^ak = 6 keV and Fx = 1.9. Unobscured narrow-line and lineless AGNS have ~100 times lower accretion rates than broad-line 
AGNs, as well as significantly cooler and somewhat weaker accretion disks. 



tively inefficient accretion flows (RIAF) expected at low 
accretion rates. The appearance of a RIAF at inner radii 
miglit produce more X-ray emission, as we discuss in §4.1 
below, but this emission probably has a similar power- 
law slope as the X-ray corona present in broad-line AGNs 
with high accretion rates. This is unsurprising, since 
both the RIAF and the corona are thought to be ionized 
plasmas with X-ray emission from inverse Compton scat- 
tering and/or bremsstrahlung. We can conclude that the 
onset of a RIAF in unobscured narrow-line and lineless 
AGNs with accretion rates of 10"'' < Lint /LEdd < 10~^ 
do not cause harder X-ray power-law slopes. 

4.1. Physics of the Accretion Disk 

As accretion rate increases from lineless and narrow- 
line to broad-line AGNs, the disk temperature signifi- 
cantly increases and its brightness with respect to the 
X-rays marginally increases. An increase in temperature 
with accretion rate is expected for a thin accretion disk, 
which has T^ax oc m^/'* (jShakura fc Sunvaevlll973D . We 
additionally discuss below how the onset of a radiatively 
inefficient accretion fiow could also cause apparent cooler 
disk emission. Both an increase in temperature and in 
Ldisk/Lx with accretion rates would contribute to the 
observed increase of aox (the ratio of rest-fra me UV to 
X-ray emission) w ith accretion rate (IKelly et al. 2008j 
Young et al.]|2010D . In our previous work (jTrump et al.l 
2009cl) . we suggested that the increase of aox with ac- 
cretion rate was due only to the disk luminosity decreas- 
ing with respect to the corona luminosity. While this is 
partly correct, the correlation is also caused by increasing 
disk temperatures at higher accretion rates. 

AGNs with Lint/ LEdd ^ 0.01 are predicted to have 
radiatively inefficie nt accretion flows (RIAFs) near the 
central black hole (iBegelman. Blandford fcReei ITOSl 
Naravan et allllQQa mianll2007lHNaravan fc McClintockl 
20081) . At such accretion rates, we can define a trun- 
cation radius Rt where the coUisional cooling time is 



comparable to the accretion time. Beyond Rt, accretion 
will remain in a standard geometrically thin and opti- 
cally thick disk with a the rmal blackbody spectrum (e.g., 
IShakura fc Sunvaevll 19731) . However within Rt, there are 
too few collisions to couple the ions and electrons and 
the gas becomes a two-temperature plasma. The elec- 
trons are cooled by bremsstrahlung, synchrotron, and 
Compton up-scattering, while the ions remain at the 
virial temperature. This means the flow is geometrically 
thick and optically thin. The introduction of a trunca- 
tion radius changes the Rin = 6Rg assumption for the 
accretion disk model, since by definition Rin > Rt- The 
peak energy of the best-fit accretion disk model is not 
very sensitive to the choice of Rin, although larger inner 
radii change the shape of the model with additional red 
emission. At accretion rates Ljnt/ LEdd ^, 10~ ^, as in 
our sample, Rt ^ 80i?g (jYuan fc Naravanll200l) . Using 
Tin = 80 in the accretion disk model fitting in Section 
3.1 doesn't change the best-flt values of Epeak, although 
it does result in slightly better fits. 

The marginal {2.1a significant) increase of Ldisk/Lx 
with Lint/ LEdd might also be caused by the onset of 
the RIAF. As Rt expands outwards, the disk emission 
decreases and the RIAF emission increases. The RIAF 
hot plasma emission is mostly X-ray bremsstrahlung and 
Compton up-scattering (like the corona), with an ad- 
ditional IR synchrotron component (which we discuss 
in Sections 4.3). As accretion rate drops and Rt in- 
creases, the rise of the RIAF X-ray emission compared 
to the optical/UV disk emission is seen as a decrease of 
Ldisk /Lx ■ Indeed, local low-luminosity AGNs have even 
lower accretion rates and larger Rt, with consequently 
lower Ldisk / Lx ratios and cooler optical thin-disk emis- 
sion (Ho"2008). 

The transition to an inner RIAF also causes the dis- 
appea rance of b road emission lines at Lint/ LEdd ^ 
0.01. INicastro I (l2000l) was the first to elegantly show 
that the broad emission lines are only present above 
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a critical accretion rate. However iNicastro"! ()200Cl[ ) as- 
su med that the innermost po ssible orbit was given by 
the ^ Shakura fc Sunvaevl (|1973[ ) thin-disk model, Vcrit — 
S.lQRg. Here we follow their basic derivation, with the 
key difference that we use the RIAF transition radius 
as the innermost orbit for the presence of a broad-line 
region. 

There is evidence that the b road emission line re- 
gion is part of a disk wind fe.g. lEmmering et al.|[T992l: 
Murrav fc ChiaMllTOQl 111^1200(11 : lElitzur fc Shlosmanl 



2006( ). The positions of individual broad emission lines 



are stratified an d set by the ionizing luminosity of the 
conti nuum (e.g. iPeterson fc Bent j l2006t iDennev et al.l 
I2009D . The base of the wind itself, however, is set by 
the radius at which the rad iation pressur e equals the gas 
pressure, defined bv .Shakura fc Sunvaevl (1973) as: 



'^wind 



(1 



-0.5 \16/21 
wind J 



15.2(aM)2/2i - 



16/21 



(8) 



with r^,nd in units of R/{6Rg) = R/{6GM/c^), M in 
units of Mbh/Mq, a is the viscosity parameter, and r/ 
is the accretion efficiency. While RIAFs are expected to 
have strong outflows (see Section 4.2), the RIAF region is 
a high-temperature ionized plasma and so any associated 
disk wind would not emit broad emission lines in the 
UV/optical. Thus the RIAF truncation radius sets the 
innermost possible radius for the existence of a broad- line 
region. Assuming that m ~ Lint/LEdd and rearranging 
Equation 8 with r^jind > Rt, ct ~ 0.1, and 77 ~ 0.1, this 
sets the minimum specific accretion rate for a broad line 
region as: 



m > 0.0U{Rt/80Rg)M'^ 



-1/8 



(9) 



with Ms = Mbh/{W^Mq). We leave Rt as a free param- 
eter since it is poorly constrained, although the best-fit 
RIAF model s for Ljnt/LEdd -^10"'^ -10"^ AGNs suggest 
Rt ~ 80Rg ()Yuan fc Naravanll2004[ ). As an AGN drops 
below this minimum accretion rate, its broad lines disap- 
pear and only narrow lines (or no lines) are observed, as 
seen in the transition at \og{Lint/ LEdd) ~ —2 transition 
i n Figures El and El 

lElitzur fc Hoi (|2009l ) also predict that the disk wind as- 
sociated with the BLR will disappear below an accretion 
rate at which the outfiowing velocity drops below the ran- 
dom velocity of the disk. Elitzur fc Ho (2009) measure 
a BLR-disappearance accretion rate of \og{L / LEdd) < 
C + B\og{Lboi) from the low-luminosity local AGNs of 
IHoI pOfM l. with ^ = -0.5 and C = 14.4. I n our sample 
(as we ll as those of iKoUmeier et al.l 120061 : iTrump et al.l 
|2009bf ). the BLR disappears at \og{L/LEdd) < 0.01. For 
a typical bolometric lum inosity of Lj^t ^ 10"*^^ '^ erg s~^ 
(also appropriate for the iKollmeier et al1l2006L sample), 
and assuming the same /3 = —0.5, this instead corre- 
sponds to C = 20.3: a remarkable difference of 6 orders 
of mag nitude. It is unlikely that the bolometric correc- 
tions of lHol (|2009| 1 are incorrect by 6 orders of magnitude, 
and so we must conclude that the Elitzur fc Ho (2009) 
model does not describe the disappearance of the BLR 
for high luminosity AGNs. Instead a disk-wind model 
following iNicastrol ()2000D best describes the BLR disap- 
pearance as the radius of wind generation region moves 
within the inner RIAF region. 



It must be noted that while disk wind models have had 
success in describin g highly ionized emission and absorp- 
tion h ues in the UV ( Proga et al.ll2000l : iProga fc Kallmanl 
l2004f ). they have not been applied to optical emission 
lines. The Ha broad emission line almost certainly forms 
in a higher density, lower ionization region than the 
Civ and Mgii broad emission lines. In addition there 
is evidence that the dynamics of the H/3 broad emis- 
sion line are wildly variable, with reverberation mapping 
indicating infalling, virialized, and outfiowing H/3 emis- 
sion regions in three AGNs (Denney et al. 2009;). While 
we do find that broad H/3 tends to be present only for 
Lint/LEdd ^ 0.01 and so fits in the wind/RIAF frame- 
work, we do not study Ha and canno t say if th is line is 
described by the same physics. Indeed. IHoI(|2009D present 
several AGNs with broad Ha emission and L^t /L Edd < 
10~^. This suggests that broad Ha emission may have 
its origin outside the disk wind, although it is important 
to note that the accretion rates of lHol ()2009l ) rely on bolo- 
metric corrections to monochromatic luminosities and so 
may suffer from significant systematic uncertainties. 

4.2. Accretion Rate and Outflows 

The gas in a RIAF is not gravitationally bound to 
the supermassive black hole because the ions are not 
losing energy through radiation. As a result, AGNs 
with RIAFs are predicted to hav e strong radio outflows 
("Nar avan et al.lll995l:lMeieHl200lD . The coupling between 
a RIAF and a strong radio outflow h as been confirmed 
by ob servations of black hole binaries ([Fender fc Bellonl 
|2004[ ). and it is poss ible to translate these observations 
to AGN scales (e.g. iMaccarone et al.l 120031 ). In Figure 
[71 we show the AGNs of our sample with the ratio of 
radio luminosity to disk luminosity. Note that since 
the radio emission is coincident with the X-ray point 
source we assume that it originates from the AGN, but 
we cannot strictly rule out other sources of radio emis- 
sion (e.g., from star formation). The Lint/ LEdd < 10~^ 
AGNs which are expected to have RIAFs tend to have 
higher ratios of radio to disk (optical/UV) luminosity. 
The mean L^isk/ Lradio for rapidly accreting broad-line 
AGNs is a factor of ten lower than the mean L^isk / L^adio 
for narrow-line and lineless AGNs, and since the scatter 
in each sample is about '-^^0.5 dex this translates to a 
highly significant difference (14.9a). 

The large scatter in the Ldisk/ Lradio ratio at both high 
and low accretion rates is likely because the radio power 
is additionally dependent on properties like black hole 
spin and orientation. But the highly significant increase 
in Ldisk/ Lradio for fow accretion rate AGNs suggests 
that Lint/LEdd < 10"^ AGNs with RIAFs generally 
have r elatively brighter radio emission. iMelendez et al.l 
(|2010f l noticed a similar trend of increasing radio lu- 
minosity with decreasing accretion rate , using [Oivjas 
a proxy for intrinsic luminosity (e.g., IMelendez et al.l 
120081: iDiamond-Stanic et al.ll2009l) . Many nearby radio 
galaxies are also measured to have low accretion rates 
and may even have their optical/UV emission domi- 
nated by synchrotron e mission rather than a thermal disk 
(jChiaberge et al.lll999D . 

In general, the radiation and disk winds of AGNs 
are thought to cause feedback on g alaxy sca l es by 
quenching star formation (e.g., iHopkins et al.l 120061 : 
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Fig. 7. — Accretion rate with a measure of radio brightness: the 
ratio of radio luminosity to disk luminosity for the 118 unobscured 
{Nh < 10^-^ cm~^) AGNs in our sample. Broad-hne AGNs are 
shown by black crosses, narrow-line AGNs by blue diamonds, and 
lineless AGNs by red squares. Narrow-line and lineless AGNs, at 
lower accretion rates than broad-line AGNs, tend to be more radio 
luminous compared to their accretion disk luminosity. 



ins fc Elvisll2010D . while radio jets are thought to 
cause larger-scale fe edback which can h eat the cores of 
galaxy clusters (e.g., iFabian et al.l 20021) and is observe d 
as extended emission line regions (iFu fc Stocktonl[2009l) . 
The fact that RIAFs tend to have stronger radio emission 
suggests that weakly accreting AGNs may remain im- 
portant for large-scale radio-mode feedback despite their 
optical/UV and X-ray luminosities. This suggests that 
heating cluster cores may not require bright q uasars, but 
can b e accomplishe d by fa int AGNs (see also lHart et al.l 
120091 ). lAllen et al.l (|2006D similarly found that several 
nearby weakly accreting AGNs had most of their Bondi 
accretion rates converted to radio outflows. 

4.3. Accretion Rate and the IR "Torus" 

A clumpy dust "torus" emits a u nique power- law signa - 
ture in the mid-IR from --1-lO^m (jisfenkova et al.ll2008[ ). 
This was first noticed observationally as a di stinct AGN 
locus in Spitzer/IRAC color-colo r space (iLacv et alJ 
[200l IStern et al.l l2005l ). although IDonlev et al.l (I2007f ) 
show that power-law selection is the most effective way 
to select AGN in the mid-IR. We compute the IR power- 
law slope in our AGNs from the host-subtracted ob- 
served IRAC photometry within the rest-frame wave- 
length range 1 < A < lO/im, shown with accretion rate 
in Figure El Type 1 AGNs typically have ajR < 0.5 
(0TR < — 0.5 in terms of the ~ form used by 
IDonlev et a l. 2007 ), matching the predictions of clumpy 
dust models (jNenkova et al.l 120081) . About 10% of Type 
1 AGNs are "hot-dust-poor" and do not satisfy the am 
selection criteriorP^. about half of the narrow-line AGNs 
and only one of the lineless A GNs lack this torus signa- 
ture. iCardamone et all (|2008[ ) similarly found that many 
X-ray AGNs did not have a mid-IR power-law, although 
they did not track it with accretion rate. In our sam- 
ple, the rapidly accreting broad-line AGNs have a mean 
ctiR = 0.38 ± 0.35, while the weakly accreting narrow- 
line and lineless AGNs have a mean am = 1.26 ± 0.72, 

For more details on this population, see lHao et al.l l|2010l ). 
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IR power-low slope a|pj (L,,, ~ u") 

Fig. 8. — Accretion rate with the power-law slope of the 1 < 
A < 10/xm IR emission for the 118 unobscured {Nh < 10^^ cm~^) 
AGNs. As in previous figures, black crosses are broad-line AGNs, 
blue diamonds are narrow-line AGNs, and red squares are lineless 
AGNs. We measure the slope ajji as L ~ i /" , correspond i ng to the 
slope /3 used in the power-law selection of IDonlev et al.l Ij20071 1 as 
/3 = ojn — 1. Most high accretion rate (iint/^'Bdd > 0.01) AGNs 
have IR power-law slopes corresponding to a dusty torus (aju < 
0.5). Of Lint / L^^^ < 0.01 AGNs, however, half the narrow-line 
and all but one of the lineless AGNs lack the torus signature. 

meaning that the two samples differ with high signifi- 
cance (10.7(t). 

A unified model based solely on geometrical obscura- 
tion suggests that narrow-line and lineless AGNs are ob- 
scure d by the same t orus present in broad-line AGNs 
(e.g.. lAntonuccil ll993D . Instead the low accretion rate 
AGNs {Lint/ LEdd < O.Ol) frequently lack the torus IR 
signature. In part, this may be because the torus power- 
law is simply being overwhelmed by the accretion disk 
SED at Lint/LEdd < 0.01. At low accretion rates, 
the temperature of the disk decreases, and a disk with 
Epeak = 1 eV will peak at 1.2 /zm, emitting a power-law 
ofQ;~2atl<A< lO^m. In a typical broad-line AGN, 
the IR to rus is roughly the sam e strength as the accre- 
tion disk (^Ric hards et al.ll2006l see also Figure [5]). Since 
many Lint/LEdd < 0.01 AGNs in Figure [8] have a > 2, 
they must be dominated by the accretion disk emission 
and have, at best, very little emission from the torus. 

The weaker or missing torus in many Lint/LEdd < 0.01 
AGNs can be described in a similar fashion to the van- 
ishing disk- wind BLR in Section 4.1. There is good ev- 
idence that the outer edge of the BLR c oincides with 
the inner edge of the c lumpy dust (jNetzer fc Laor. ,19931: 
iSuganuma et al.l[2006f l. Some authors additionally sug- 
gest that the BLR and the clumpy dust "torus" are two 
compone nts of the same wind driv en off the accretion 
disk fe.g.. lElitzur fc Shlosman|[2Q06l ). If the clumpy dust 
wind emerges from the disk at a similar radius to that 
calculated in Section 4.1, then we would expect the IR 
power-law signature to disappear at Lint/LEdd < 0.01, 
just as the BLR disappears. However many narrow-line 
AGNs with Lint/LEdd < 0.01 still have the negative 
IR power-law slopes, suggesting that there must be an- 
other source of mid-IR emission. Either there is a distant 
source of clumpy dust beyond the expanding RIAF, or 
there is mid-IR synchrotron emission in th e RIAF region 
at th e base of the radio jet (as observed bv lLeipski et al.l 
l2009h . 
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5. A SIMPLE MODEL FOR UNIFYING AGNS BY SPECIFIC 

ACCRETION RATE 

Figure [5] presents a simple schematic outlining the 
changes in AGNs from high {Lint/ LEdd > 0.01) to low 
{Lint/ LEdd < 0.01) accretion rate. At the top is a broad- 
line AGN with high accretion rate {Lint/ LEdd ~ 0.1). At 
these high accretion rates the gas and dust falling into 
the black hole forms a thin accretion disk and a disk 
wind originates at Rwmd ^ 250i?g. The broad emission 
lines are emitted in stratified regions along this wind 
based on the radiation pressure (which ionizes and ex- 
cites the wind material), with Rblr ^ L^-^ and high 
ionization lines (e.g., Civ) emitted from nearer radii 
than low ionization lines (e.g., Mgli) (jPeterson fc Bent3 
l2006f). At higher radii, the disk wind forms clumpy dust 
(jNenkova~et al. 200_8i)- This dusty "torus" can otascure 
the AGN along lines of sight near the disk, causing an 
obs erver to see an obscured narrow-line AGN (jAntonuccil 
[199^. 

The bottom of Figure [9] shows an AGN with low ac- 
cretion rate {Lint/ LEdd ~ 0.003), characteristic of the 
unobscured narrow-line and lineless AGNs in our sam- 
ple. The onset of a geometrically thick RIAF changes 
the picture dramatically. Because the disk wind radius 
is within the RIAF, there are no broad emission lines. 
Instead the dominant outflow is a radio jet, and AGNs 
with low accretion rates and RIAFs are typically more 
radio luminous than broad-line AGNs. The lack of a disk 
wind also means that there is not the typical clumpy dust 
"torus" seen in broad-line AGNs. However we cannot 
rule out the presence of dust completely, as clumpy dust 
may come from another source besides the disk wind and 
some Lint/ LEdd 0.01 have the IR signature of hot dust. 

6. PREDICTIONS AND FUTURE OBSERVATIONAL TESTS 

The multiwavelength data of COSMOS provide many 
diagnostic capabilities, and we have argued that de- 
creasing accretion rates lead to the onset of a RIAF 
at m < 0.01 and subsequently stronger radio jets, a 
weaker torus, and the disappearance of broad emission 
lines. The onset of a RIAF also makes several predic- 
tions testable by future observations. In addition the 
simple model in Section 5 can be more fully constrained 
by additional investigations. 

If the broad-line region is truly disappearing at m < 
0.01 then we would expect spectropolarimetry to re- 
veal reflected broad emission lines in only high accre- 
tion rate (m > 0.01) narrow-line and lineless AGNs. 



Spectropolarimetry of nearby AGNs shows a dichotomy 
based on accretion rate, although most authors place the 
change from hidden broad lines to "true" Type 2 AGN 
at m ~ 0.001 (|Tranl 120031 : 1 Wang fc Zhand [2007[ ) . Most 
likely, the difference results from the uncertain bolomet- 
ric corrections used in these previous works, compared 
to the full modeled SEDs used here. 

Mid-IR broad-band polarimetry could determine the 
cause of the negative IR power-law slopes in m < 0.01 
AGNs. If the clumpy dust "torus" is associated with the 
same wind that drives the broad line region, it should 
disappear in these objects. The mid-IR signature might 
instead be synchrotron radiation in the RIAF at the base 
of the jet, which would appea r polarized at the > 3% 
level (e.g.. iJannuzi et anil994[ ). If no polarization is de- 
tected, then we must conclude that clumpy dust exists at 
higher radii than the BLR disk-wind, beyond the RIAF 
region of m < 0.01 AGNs. 

It is very difficult to measure accretion rates of par- 
tially or fully obscured AGNs, and such objects are gen- 
erally missed by the X-ray and optical limits of this study. 
However we do make a few predictions for the accretion 
rates of various AGNs. If the torus is part of a disk-wind 
that vanishes at rh < 0.01, the n torus-obscured AGNs 
of the classical I Antonuc"dl (|1993[ ) unified model will have 
only high accretion rates (m > 0.01). Obscuration by 
cooler dust associated with host galaxy star formation, 
as predicted by the observed reds hift evolution in the 
narrow- line/broad-hne AGN ratio (jTreister et al.l 120091 : 
iTrump eFal... 2009a) , could conceivably be present at any 
accretion rate (although it may be limited by the abil- 
ity of the dusty star formation to feed the black hole, 
lBallantvngl2008[ ). We might then expect that obscured 
AGNs with a strong mid-IR torus signature should have 
m > 0.01, while AGNs obscured by the cooler dust as- 
sociated with host galaxy star formation might have a 
wider range of accretion rates. 
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Fig. 9. — A schematic model showing the changes in the accretion disk from a broad-line AGN with high accretion rate {Li„t / I^Edd ~ 0.1) 
to a narrow-hne or hneless AGN with low accretion rate {Li„t / LEdd ~ 0.003). The x axis shows the radial distance from the black hole in 
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^ The full catalog will appear as a machine-readable table in the electronic version. 

^ "BL" refers to a broad-line AGN, "NL" is a narrow-line AGN, and "OD" is a lineless or optically 

dull AGN. 

"S" means the spectrum a nd redshift are from the SDSS archive, "I" is from the COSMOS 
Magellan/IMACS campaign (jTrump et al.ll2009al 'l. and "Z" is from the zCOSMOS VLT/VIMOS 
campaign (,Lillv et al.ii2007i') . 



